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Abstract. The usefulness of study of hyperfine splitting in the hydrogen atom is 
limited on a level of 10 ppm by our knowledge of the proton structure. One way to 
go beyond 10 ppm is to study a specific difference of the hyperfine structure intervals 
8Ai'2 — Avi. Nuclear effects for are not important this difference and it is of use to 
study higher-order QED corrections. 



1 Introduction 

The hyperfine splitting of the ground state of the hydrogen atom has been for a 
while one of the most precisely known physical quantities, however, its use for 
tests of QED theory is limited by a lack of our knowledge of the proton structure. 
The theoretical uncertainty due to that is on a level of 10 ppm. To go farther 
with theory we need to eliminate the influence of the nucleus. A few ways have 
been used (see e. g. ||]): 

• to remove the proton from the hydrogen atom and to study a two-body 
system, which is like hydrogen, but without any nuclear structure, namely: 
muonium j2| or positronium 

• to compare the hyperfine structure intervals of the Is and 2s states (this 
work); 

• to measure the hyperfine splitting in muonic hydrogen and to compare it 
with the one in a normal hydrogen atom (a status report on the n — 2 
muonic hydrogen project is presented in Ref. Q]; comparison of the Is and 
2s hfs and possibility to measure Is hfs is considered in Ref. 

Recently there has been considerable progress in measurement and calculation 
of the hyperfine splitting of the ground state and the 2s 1/2 state in the hydrogen 
atom. The 2si/2 hyperfine splitting in hydrogen was determined to be |^ 

Z\i^2(H) = 177 556.785(29) kHz , (1) 

While less accurate than the classic determination of the ground state hyperfine 
splitting, the combination of Is and 2s hfs intervals 
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which is determined in the hydrogen atom Q as 

D2i{H) = 48.528(232) kHz , 

has more impUcations for tests of bound state QED because there is significantly 
less dependence on the poorly understood proton structure contributions. Specif- 
ically, the theoretical uncertainty for the ground state from the proton structure 
is about 10 kHz, while the uncertainty for the combination is estimated to be 
few Hz. 

On the theoretical front, there has been considerable progress in the calcula- 
tion of the ground state hyperfine splitting. Taken together with earlier calcula- 
tions of Dzi 0,111, which were possible because of cancellations of a number of 
large terms, one can now give quite accurate values for Ai^i and Av2, We collect 
in Tables 1 and 2 along with the hydrogen results, the known experimental and 
theoretical results for the deuterium atom and the ^He"*" ion. The helium results 



Avi{^Ile+) = 8665 649.867(10) kHz (3) 

and § 

Z\i'2(^He+) = 1083 354.969(30) kHz (4) 
lead to the most accurate value for the difference 

i:»2i(^He+) 1 189.979(71) kHz . (5) 



Table 1. Comparison of the QED part of the theory to the experiment for hydrogen 
and deuterium atoms and for the '^He^ ion. The results are presented in kHz 



Atom Experiment QED theory for D21 





D2i(exp) 


Refs.: 2s/ Is 


Old 


New 


H 


48.528(232) 


§/§ 


48.943 


48.969(2) 


H 


49.13(40) 


0/i 






D 


11.16(16) 


0/0 


11.307 


11.3132(4) 


3He+ 


1 189.979(71) 


i/i 


1 189.795 


1191.126(40) 


3He+ 


1 190.1(16) 


iii/i 







2 Theory 

We consider a hydrogen-like system with a nucleus of charge Z , mass M, spin /, 
and magnetic moment /i. The basic scale of the hyperfine splitting is then given 
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by the Fermi formula, 

E,^'-Z^a^Ryd\A'-^(^\\ (6) 
^3 " ^iB 21 \m + Mj ^ ' 

Here we take the fine structure constant a derived from g-2 value of electron 
a^^ ~ 137.035 999 58(52). In addition we use a value the Rydberg constant of 
Ryd = cRy^ 3.289 841 931 • lO^^ kHz. 

We present the hyperfine structure as a sum 

Avn — Avn{QED) + Z\i/„ (nuclear structure) . (7) 



2.1 Non- recoil limit 



First we consider the external-field approximation. For a point-like nucleus, they 
can be compactly parameterized by the equation 



Z\^.„(N-R) = ^ 



Here, with 7 — i/l — (Za)^, |^ the Breit relativistic contribution is 

1 



Bi = 



7(2,-1) + + ' 



and 



2(2(l + 7) + v/2(l + 7)) 1 ^ 17,7 ,2^ 449 4 

-D2 = — ^ tt; — r- — ^ — - — 1 H (Za) H (Za) 

(1+7)27(472-1) 8 ^ ^ 128^ ^ 



(8) 



(9) 



and the functions D^[Za) represent r loop radiative corrections. In the limit 
Za = they reduce to the power series expansion of the electron <? — 2 factor, 
and the difference is refered to a binding correction. For the ground state. 



D 



(2) 



1 



+ 7r(Za) ln(2) 



[Zaf 



;\Tl^{Za) 



V 3 

and for the excited state 



( }^ in(2) - ^ ) In(Za) + Gf'iZa) + G^^iZa) 



180 



(10) 



Df' = -+TT{Za) ( ln(2) - - ) + (Za 



1 



ln^(ZQ;) 



32 



ln(2) 



1541 
180 



In(Za) + GT{Za) + GY{Za) 



(11) 
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At present the functions G^^ have been determined numericaUy at Z = 1 and 

Gf^{Z = 1) = 16.079(15) 

and 

Gf^{Z = 2) = 15.29(9) 
while GX^ is known analytically |p^ : 



(12) 
(13) 



Gl 



VP 



15 



ln(2) 



34 
225 



+ TT(Za) 



13 Za 

24 



539 
288 



(14) 



To present results for the 2s state, we can use the results of Ref. ||lj| for 1321, 
which however include terms only up to order a{Za)'^Ep. After we recalculated 
some integrals from paper |lj] the result is 

16 



7 + — ln(2) In(Za) - 5.221233(3) + 0(7r(Za)) (15) 



and 



Gr = Gr-^ + ^ln(2) + 0(7r(Za)) 



(16) 



Continuing to the two-loop corrections, all terms known to date are state- 
independent, so we give only the ground state result | |T^ , p^ , pO| , 

l(4) 



D\^> = a(4) + 0.7718(4) 7T{Za) 



[Zaf\n{Za) 



Non-leading terms, including single powers of ln(ZQ;) and constants, are both 
state-dependent, but unknown. 

When the nucleus is not point-like, the leading correction is known as the 
Zemach correction. 



Z\t'„(Zemach) 



^Ef 



{Zam) j ^ {Ge{p)Gm{p) 



1 



(17) 



Inaccuracy arisen from uncertainties in the form factors Ge and Gm, which are 
both normalized to unity at zero momentum, and from the lack of knowledge of 
polarization effects, is large as about 10 ppm or 4 ppm respectively, but those 
leading terms are state-independent and do not contribute into the the difference 

For atoms with nuclear structure the following result was found p5| 



AD2i{Rec) = (Zaf 



9 



M 



7 

32 



ln(2) 



1 



145 7 , 

ln(2 

128 8 ^ ' 



(18) 



where gM/Znip = x = {fi/ns) (M/m) (l/ZI). It does not depend on the 
nuclear structure effects such a distribution of the nuclear charge and mag- 
netic moment. Contrary, the leading recoil term for the Ai>n (which has order 
{Za){'m/M)hi{M /m) pl| is essentially nuclear-structure dependent. 
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3 Present status of D21 theory 
3.1 Old theory and recent progress 

The Breit, Zwanziger and Sternheim corrections p^ , p^Jl5| lead to a result 

"5 177 



D21 = Ep (Zayx 



128 



a 

TT 

a 

TT 

m 
M 



(Za) 



16 



-7 + y ln(2) ) In(Za) - 5.37(6) 



9 



7 

32 



ln(2)\ 



1 

1 - - 

X 



/ 145 
V128 



ln(2)) 



(19) 



Some progress was achieved before we started our work. In particular, we 
need to mention two results: 

• Integrals, used for in Ref. 10, were evaluated later by P. Mohr^ with higher 
accuracy and the constant was found to be -5.2212 instead of -5.37(6). The 
theoretical prediction based on Eq. ( |l9| ) but with a corrected value of the 
contstant is Table 1 as "old theory" . 



• Some nuclear-structure- and state-dependent corrections were found 
arbitrary nS. 



for 



3.2 Our results 

The similar difference has been under investigation also for the Lamb shift and a 
number of useful auxiliary expressions have been found for calculating the state 
dependent corrections to the Lamb shift |Q . 

Let us mention that an improvement in the accuracy and new result on higher 
n hfs can be expected with progress in optical measurements and we present here 
a progress also for higher n, defining Dni — n? Avn — Avi. 

• We have reproduced the logarithmic part of the self energy contribution and 
found for arbitrary ns 



ADni = - {ZaY Ep In(Za) 

TT 



- ln(n) 



^ Unpublished. The result is quoted accordingly to Ref. 



2n2 



(20) 
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The calculation is based on a result in Ref. for the singe logarithmic 
correction due to the one-loop self energy and one-loop vacuum polarization. 



We have reproduced the vacuum polarization contribution and found that 
for arbitrary ns 



TT V 15 



TL — 1 

2 ( - ln(n) H h VX") " '0(1) 

11 



2vP- 



(21) 



Integrals used by Zwanziger jlj] have been recalculated and the constant 
was found to be -5.221233(3). 

We also found two higher-order logarithmic corrections 



and 



ZiD„i = {ZaY Ep \n{Za) 



n — 1 

2 ( -ln(n) H h ■0("-) --0(1) 



- 1 
2^2 



(22) 



OL m , „ -.2 ^ ' ^ N /16 



Z\Aa = -77 {ZaY Ef In(Za) — 



n - I 

2 ( - ln(n) + V(") + V'(l) 



- 1 



2n2 



(23) 



We found two higher-order non-logarithmic corrections 

139 



16 



41n(2) 



n - 1 

-ln(n)H |-0(n) - V'(l) 



ln(2) - H 



ip{n + 1) - ■0(2) - ln(n) 



(n- l)(n + 9) 



(24) 



and 



5 

24 



Ti — 1 

-ln(n) H h'0(n) - ■0(1) 
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+ - X tp{n + 1) - tp{2) - ln(n) - 



(n-l)(n + 9) 
4n2 



1 



(25) 



• We have also found a term proportional to the magnetic radius. To the best 

of our knowledge that is the first contribution, which is proportional to the 
magnetic radius and on the level of the experimental accuracy. The complete 
nuclear-structure correction is 



4 Present status 

To calculate the corrections presented in the previous sections, wo liavc used 
an effective non-relativistic theory. In particular we have studied two kinds of 
terms. One is a result of the second order perturbation theory with two S{r)- 
like potentials, evaluated in the leading non-relativistic approximation, while 
the other is due to a more accurate calculation of a single delta-like poten- 
tial. Both kinds contribute into the state-independent leading logarithmic cor- 
rections {a'^{Za)'^\n^{Za), a{Zaf{m/M)\T?{Za), and a{Z af \ti{Z a)) and to 
next-to-leading state-dependent terms {a'^{Za)'^ ln{Za), a{Za)'^{m/M)ln{Za), 
and a{Za)^). The crucial question is if we found all corrections in these orders. 
Rederiving a leading term in order a{Za)^ \n{Za) within our technics, we can 
easily incorporate the anomalous magnetic moment of the electron in the cal- 
culation and restore the nuclear mass dependence. Since we reproduce the well- 
known result for the a{Za)^ \n{Za) term, we consider that as a confirmation of 
two other logarithmic corrections found by us. In the case of a(Za)^ it might 
be a contribution of an effective operator, proportional to {A/m)S{r). That can 
give no logarithmic corrections, but leads to a state-dependent constant. We are 
now studying this possibility. 

Summarizing all corrections, the final QED result is found to be: 



^^Nuci ^ _^zaf ^(n + 1) - ^(2) - ln(n) - 



(n-l)(n + 9) 
4n2 



X Ziz^i(Zemach -|- polarizability) -|- - {Za)"^ ^^{n) — ?/^(l) — ln(n) 




(26) 
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a 

+ - 

TT 



_1 + Aln(2) 
10 15 ^ ' 




ln(2) 



ln(2)j In(Za) 
ln(2) ) \n{Zo 



139 5 

4 n(2 H 

16 ^ 24 



145 
128 



ln(2) 



- ln(2) 



13 

T 



ln(2) 



ln(2) 



3 
16 



(27) 



Numerical results (in kHz) for hydrogen and deuterium atoms and the helium- 
3 ion arc collected in Table 2. One can see that the new corrections essentially 
shift the theoretical predictions. A comparison of the QED predictions (in kHz) 
against the experiments is summarized in Table 1. We take the values of the 
fundamental constants (like e. g. the fine structure constant a) from the recent 
adjustment (see Ref. [p5|). 



Table 2. QED contributions to the D21 in hydrogen, deuterium and helium-3 ion. The 
results are presented in kHz 



Contribution 


H 


D 


3He+ 


[ZafEp 


47.2275 


10.8835 


1 152.9723 


+ a{ZafEF (SE) 


1.9360 


0.4461 


37.4412 


+ a{ZafEF (VP) 


-0.0580 


-0.0134 


-1.4148 


+ {Zaff^EF 


-0.1629 


-0.0094 


0.7966 


+ a^{ZafEF 


0.0033(16) 


0.0008(4) 


0.070(35) 


+ a{ZaffjEF 


-0.0031(15) 


-0.0004(2) 


-0.022(11) 


+ a{ZafEF (SE) 


0.0282 


0.0065 


1.3794 


+ a{ZafEF (VP) 


-0.0019 


-0.0005 


-0.0967 



An important point is that the difference is sensitive to 4th order corrections 
and so is competitive with the muonium hfs as a test of the QED. The difference 
between the QED part of the theory and the experiment is an indication of 
higher-order corrections due to the QED and the nuclear structure, which have 
to be studied in detail. In particular, we have to mention that while we expect 
that we have a complete result on logarithmic corrections and on the vacuum- 
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polarization part of the a{Za)^ term we anticipate more contributions in the 
order a{Za)^ due to the self-energy. A complete study of this term offers a 
possibility to determine the magnetic radius of the proton, deuteron and helion- 
3. 
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